Strong field effects in laser pulses: the Wigner formalism. 
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We investigate strong field vacuum effects using a phase space approach based on the Wigner for- 
malism. We calculate the Wigner function in a strong null-field background exactly, using lightfront 
field theory. The Wigner function exhibits the distinct features of strong field QED: in particular 
we identify the effective mass in a laser pulse, compare it to the well-known mass shift in a periodic 
plane wave and identify signals of multi-photon absorption and emission. Finally, we show how to 
extend our results to describe vacuum pair production in colliding laser pulses. 
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I. INTRODUCTION. 

Vacuum phenomena generated by strong external 
fields are currently a topic of considerable interest, due to 
the prospects provided by infrastructure projects such as 
the European XFEL Q and the ELI facility @,||]. The 
latter is expected to produce laser pulses with an inten- 
sity of 10 24 W/cm 2 and peak electric field strengths of 
10 14 V/cm. These fields can produce strongly relativis- 
tic ionisation effects, polarise the vacuum and so gener- 
ate nontrivial vacuum refractive indices , be used 
to search for new light particles Q, and, perhaps, pull 
electron-positron pairs out of the vacuum [8[. 

Such strong field physics has been investigated using 
scattering calculatio ns p , exact solutions semiclas- 
sical techniques [TT| - [T^ ]7 Monte Carlo simulations [l5[ 
and quantum kinetic equations [l6|-[l8|. A framework 
which includes quantum kinetic theory is provided by the 
Wigner function, the quantum analogue of a probability 
distribution. It is a well-established tool in quantum op- 
tics [l9|], and has been applied to transport phe nomena 
in plasmas and condensed matter systems j20ll2l| , and to 
pair pr oduction in both strong electric and colour fields 
[22H24j . (For connections with other approaches to pair 
production, see [25|,[26[.) The Wigner formalism has the 
advantage of offering a phase-space description of quan- 
tum systems with direct access to real time processes. 
For a thorough introduction see [27I . 

One of the challenges of this approach lies in its in- 
terpretation: quantum effects can generate negative val- 
ues of the Wigner function which have no probabilis- 
tic interpretation. Nevertheless, the Wigner function is 
an observable, and indeed has been measured in vari- 
ous (non-rclativistic) optical and atomic setups [28| - |30| . 
Consequently, it is of great importance to obtain a thor- 
ough analytical understanding of the Wigner function. In 
this paper we therefore calculate the exact (relativistic) 
electron-positron Wigner function, in the presence of a 
strong background field, and give a complete description 
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of its properties. Our chosen backgrounds, modeling high 
intensity lasers, arc null-fields, i.e. plane waves. This au- 
tomatically includes effects due to strong magnetic, as 
well as electric, fields. Although we neglect transverse 
size effects, we take explicit account of effects due to finite 
pulse duration by considering finite longitudinal extent of 
the fields. 

The most well known effect of a plane wave on an elec- 
tron is the generation of an intensity dependent effective 
mass [311 ]. This is only well understood for very long 
laser pulses modelled by periodic plane waves (i.e. those 
of infinite duration), where rapid changes in the parti- 
cle's energy and momentum, due to oscillations of the 
fields, average out to a mass shift. Modern laser pulses 
are, however, only a few cycles long. In order to exam- 
ine the impact of finite beam geometry, we will compare 
results for periodic plane waves with those for more re- 
alistic pulse profiles. Using the Wigner function, we will 
identify the effective mass for an arbitrary profile. 

This paper is organised as follows: wc begin in Sect. |TT] 
by briefly reviewing the gauge invariant Wigner func- 
tion in both covariant and equal time formulations. Wc 
then recap the behaviour of charges in a plane wave, in 
Sect. IIIII In Sect. IIV1 we construct the Wigner function 
in a plane wave background exactly, and identify the ef- 
fective mass in a pulse. We examine in detail the Wigner 
function for both oscillatory and pulsed fields, identifying 
effects due to the effective mass, multi-photon contribu- 
tions and quantum corrections. Finally, we add a second 
laser pulse and describe how our results may be extended 
to include pair production. 



II. THE WIGNER FUNCTION. 



The electron-positron Wigner function in an external 
electromagnetic field F^ v is constructed from the field 



density in the vacuum state | ) : 



f)]|0) 



=:Z4(z+f,z-§), (1) 



where x and y are the 'centre of mass' and 'relative' co- 
ordinates, respectively. The covariant Wigner function 
itself, W4, is obtained from U4 by Fourier transforming 
with respect to the relative co-ordinate, trading y^ for a 
momentum whose physical interpretation will be given 
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directly: 



W 4 (x,p) = / d 4 y e }p-y/K e ieSdz.A(z)/h Ui 



(2) 



The Wilson line in ^ makes W4 gauge invariant and is 
evaluated along the straight line path from x — y/2 to 
x + y/2 [20I ] . Expanding in powers of ft, one can show 
that the classical part of W4 is always constrained to obey 
p 2 — m = 0, and hence p^ may be identified with the 
kinematic momentum. As we will calculate exactly in h, 
we set it equal to one from here on. The Wigner function 
obeys a manifestly gauge invariant version of the Dirac 
equation 27j, which we will return to in Sect. |Vj Conse- 
quently, there are at least two approaches to construct- 
ing Wigner functions. First, by solving this differential 
equation, which is a manifestly gauge invariant approach. 
Alternatively, one can fix the gauge, quantise the theory 
and construct the expectation value in ((4]) directly. We 
give here a simple example of the latter method: the free 
Wigner function is easily found by quantising the free 
fcrmion theory and calculating the commutator in the 
free vacuum state. One obtains 



W 4 (x,p) h . cc = 2tt S(p 2 -m 2 ) 
= W 4 (x, p) f Iee [p 



(3) 



so that the free Wigner function is spatially homogeneous 
and tells us that free particles are on-shell. The final 
equality indicates a separation of W4 into spin and scalar 
contributions - in fact, W4 is the Wigner function for 
scalar particles, as is easily found by calculating the free 
scalar density 

U 4 (x + §,z-§):= (0 \{</>{x + f), 4>\x - f )}| 0) , (4) 

and performing the analogue of the transformation ([2]). 

One often also considers the equal time Wigner func- 
tion W3. This is defined as the energy average of W4, i.e. 
it is obtained by integrating out po. For the free theory 
one therefore has: 



W 3 (t;x,p) freo := / ^ W 4 (a;,p) frec 



y/p 2 



(5) 

At this point we introduce lightfront co-ordinates = 
x° ± x 3 , x 1 - = {x x ,x 2 } which will be useful below. In 
lightfront field theory (for reviews see [H, [H| and refer- 
ences therein) one quantises at equal lightfront time x + . 
One simplification of this approach is that only half of a 
Dirac spinor's components represent dynamical degrees 
of freedom [32|, Introducing the orthogonal projec- 
tion operators A* = 4 -7 :F 7 ± , the dynamical fermion field 
is = A + ?/>. On the other hand, tp(-) = A~-0 is a 
constrained field which is determined by ip(+)- I n this ap- 
proach we therefore construct the Wigner function from 
the commutator of the dynamical field "0(+) m the light- 
front vacuum. The density becomes 



and the Wigner function is again found via ([2]). The 
spin structure in this picture leads to a very simple re- 
lation between the spinor and scalar Wigner functions: 
one finds that the free Wigner function is given by the 
scalar result multiplied by the physical projector, so 



W 4 (x,p) frcc = 2p_A + W4(x,p){ le 



(7) 



where, of course, the free scalar Wigner function matches 
that found in the usual approach, see ©. As we will see, 
relation ([7]) also extends to the interacting theory. For 
future use we note that, using lightfront co-ordinates, the 
free scalar Wigner function may be written 



W 4 (x,p) ilee 



2\p 



Slp + - P ^-), (8) 



4p„ 



which simply exhibits the mass-shell condition in light- 
front co-ordinates. It follows that the equal lightfront 
time Wigner function in the free fermion theory is (writ- 
ing x = {x-,x ± }, p = {p_,p ± }) 



2n 4|p„| 



W 3 + 0r + ;x, p) fre 

/ Z7T i\p 

(9) 

where the '•' simply illustrates the separation between 
spin and scalar contributions. 



III. PLANE WAVES AND THE MASS SHIFT. 

Wc now review the effects of a plane wave background 
on particle motion and QED scattering events. The field 
strength of a general plane wave is represented by 



(10) 



)]|0) 



(6) 



where k 2 — 0, j is summed over the two transverse direc- 
tions and a dot is differentiation with respect to k.x (for 
convenience). The polarisation vectors obey a l .k = 0, 
a 1 .a J = — (rr? a 2 1 e 2 )b % \ which defines an invariant, di- 
mensionless amplitude a. Since fc M is a lightlike direction 
we take fc M = w(l,0,0, 1) so that k.x = lox + is propor- 
tional to lightfront time. 

Consider first the Lorentz equation for the classical 
orbit of a charge in this background. This is exactly 
solvable: k.p is conserved and hence the charge's proper 
time is proportional to lightfront time x + . In long pulses 
modelled by periodic fields, i.e. plane waves of infinite 
extent, the average, or 'quasi' momentum of the charge 
over a laser cycle is of particular importance: it is on- 
shell, but for a shifted mass m*, i.e. (p) 2 = m 2 . For 
periodic, infinite plane waves - IPWs from here on - the 
mass shift is m 2 = m 2 (l + aQ), where we have introduced 
the common measure of laser intensity Oq, which equals 
a for circular polarisation and a 2 /2 for linear polarisa- 
tion [13 ■ Currently, a - O(10), while for ELI-strength 
lasers a ~ 10 3 - 10 4 . 

The effective mass is only well-understood in an IPW, 
where the relevant average is taken over a laser cycle. 
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FIG. 1. Emission spectrum for photons scattered from an 
electron-laser collision, as a function of scattered photon fre- 
quency v' — w'/m. Top: n = 1 to 4 laser photon contribu- 
tions (harmonics) at an intensity ag = 25, compared to the 
ordinary Compton scattering spectrum (n = 1 at ao = 0). 
Bottom: the full spectrum for oq — 25 and a% = 1/2, summed 
over all photon numbers. 



The IPW model was assumed in the analysis of the SLAC 
E144 pair production experiment [HJ], in which the laser 
pulse contained around 10 3 cycles. To model modern 
femtosecond duration pulses, it is necessary to look at 
pulse profiles with finite extent in k.x. Periodicity is then 
lost, and thus it is unclear what average momentum is 
important, nor if any mass-shift remains. We will address 
this shortly. Suppose for a moment that some average 
( ) over lightfront time can be identified, be it the cycle 
average or otherwise. Then, from the classical orbit, one 
finds 



(p) 2 = m 2 (l + 



(11) 



where var 2 [/] = (f 2 ) — (f) 2 . Thus, the mass shift in 
periodic fields is given by the variance of the integrated 
field strength over a laser cycle. It does not seem that 
this statement, which will be useful below, is well known, 
although we believe the first appearance of the variance 
expression is in [39| . 

Going to QED, one finds that scattering amplitudes 
have kinematic support on the conservation of quasi- 
momcntum: hence, for example, electron-positron pair 
production stimulated by incoming probe photons is de- 
scribed in terms of the production of mass m* pairs [35F - 
l37j . In addition to mass-shift effects, QED scattering 
processes in periodic plane waves exhibit effects due to 
multi-photon absorption from the background. In non- 
linear Compton scattering, i.e. the emission of a photon 
from an electron in a laser field, the mass shift is respon- 



sible for a red-shift of the kinematic Compton edge, while 
multi-photon effects allow for higher harmonics, that is 
scattered photons with energy higher than the edge value 
Q. The full scattering process is a sum over all multi- 
photon contributions. This is illustrated in Fig. [TJ (Sec 
40 42] for corrections due to finite size effects.) This 
completes our discussion of the most well-known effects 
in periodic plane waves. We now turn to a closer inves- 
tigation of both periodic and pulsed laser field effects, 
using the Wigner function. 



IV. AN EXACT WIGNER FUNCTION 

In this section we calculate the Wigner function in a 
plane wave field exactly. The Klein-Gordon and Dirac 
equations are exactly solvable in this background [43| , at 
least when one works in lightfront gauge k.A ~ A + = 0, 
when the gauge potential for our F^ v becomes 



A^k.x) = fj(k.x)a J 



(12) 



In this setup, we have a correspondence with, e.g., [23J, 
which takes to be a spatially homogeneous, time de- 
pendent electric field. Instead of ordinary time t, we have 
front time x + . Instead of temporal gauge, A = 0, we 
have lightfront gauge A + = 0. The t dependent electric 
field is replaced by x + dependent electric and magnetic 
fields, and all fields are homogeneous in their correspond- 
ing 3-spaccs. The real difference is that the former field 
produces pairs, the latter does not. We will include a 
second plane wave, with which we can generate pairs, in 
a later section. 

To proceed, one quantises at equal lightfront time and 
calculates the field density in the lightfront vacuum state. 
Just as in the free case (J7J), one finds that the electron- 
positron Wigner function is given by the a scalar contri- 
bution multiplied by the background-independent spin 
contribution 2p_A + . For this reason, we restrict the re- 
maining discussion to the interesting scalar part (which 
is again the Wigner function for scalar fields) . One finds 
that the scalar U4 is 



d 4 q 
(2tt) 3 



x exp 



2k.q 



6(q 2 - m 2 )e- iq - y 

k.(x+y/2) (13) 

ds [2eA(s).q-e 2 A 2 (s)} , 



k.(x-y/2) 



where we recognise the usual Volkov phase factor [43J. 
It is possible to perform the y~ and y^ integrals in © 
exactly, to obtain the compact, and exact, expression 



Wa{x, P ) 



1 



4|p_ 



dy + exp 



iy \p+ 



pj + M 2 
Ap_ 



(14) 
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where M 2 (k.x, k.y) := m 2 (l + a 2 var 2 [/j]), with the aver- 
age taken over the separation of the fields, i.e. 



k.(x+y/2) 



ds 



k.(x-y/2) 



m 

k.y 



(15) 



When the fields turn off, M — > m and the y + integration 
correctly recovers the delta function ((5J), imposing the 
free mass-shell condition. Hence (|14[) describes a defor- 
mation of this condition in which the free mass is replaced 
with an effective mass M of the IPW form (TTTT) : hence, 
our Wigner function gives us the explicit form of the effec- 
tive mass for an arbitrary pulse profile. This is corrobo- 
rated by [44[ , where precisely this M appears in the gauge 
invariant part of the electron propagator. (See [451 ] for 
related non-relativistic results.) Let us compare the ef- 
fective mass in a circularly polarised IPW (f\ = sin (k.x), 
/a = cos (k.x)) with that in an n-cycle pulse, where [46| 



fi(k.x) = sin 4 (|^) sin(k.x) 
f 2 (k.x) = sin 4 (|f )cos(fe.a;) 



(16) 



for < k.x < 27m and zero otherwise. A pulse with 
n = 8, for example, corresponds to 32 fs lab duration at 
optical frequency. The effective mass in the IPW is 



M c 2 irc /m 2 = 1 + a 2 - a 2 sinc 2 ( 



k.y ' 

2 > 



(17) 



which is independent of k.x (linear polarisation intro- 
duces periodicity in k.x). We do not write down the 
lengthy expression for the pulsed M explicitly. The ef- 
fective mass in both the IPW and the pulse (fTB"|) are plot- 
ted in Fig. The IPW result is shown in the top panel, 
with the interesting behaviour in the k.y direction: at 
k.y = 0, the effective mass is the free mass. As k.y ap- 
proaches 2-7T, i.e. a front time duration corresponding to 
one laser cycle, M 2 increases quickly toward m 2 (l + a 2 ), 
and converges toward this value as k.y —> oo. This is 
of course the IPW m 2 associated with the cycle average: 
hence, as we know, a particle seeing infinitely many laser 
cycles behaves as if its rest mass were m*. 

In the bottom panel, we see that the effective mass in 
a pulse behaves quite differently. For large k.x, M —> to,, 
as one sits far outside the pulse. As the front time sepa- 
ration k.y increases from zero, M rises inside the pulse, 
but then falls back to zero for increasing k.y. This re- 
flects the statement that any massive particle in the field 
will leave the pulse after a finite time (front or instant), 
so that eventually the effective mass reduces to the free 
mass. In particular, we see rather cleanly that scattering 
processes in a pulse, which are asymptotic in k.y, must 
be sensitive to the free mass. 




FIG. 2. The mass shift M 2 /m 2 - 1 in units of a 2 , i.e. var 2 [/,]. 
Top: in a circular polarised IPW the mass shift tends to 1, 
i.e. M 2 — > m 2 , as k.y increases. Bottom: in the 8-cycle pulse 
(fl6)) . M reduces to the free mass m for large k.x or k.y. 



under the integral of (fl~4"|) and take e — > 0. (We have con- 
firmed that our results are unchanged using a different 
regulator.) Like the free solution, our W4 is spatially 
homogeneous, as follows from (|17l) being independent of 
k.x. W4 is plotted in Fig.[3]as a function of p_, with other 
variables given in the caption^. First we demonstrate the 
convergence of IF4 for different values of e: convergence is 
very rapid for nearly all p_ values as e decreases. Slower 
convergence is seen around several sharp edges in W4, 
which will be explained very shortly. The final result 
for W4 is represented by the solid line, which we now 
analyse in more detail. While the free Wigner function 
is strongly peaked (with a regulator understood) around 
the free mass-shell condition p 2 = to 2 , the strong peak 
in Fig. [3] sits at the shifted mass-shell condition p 2 = m 2 
of the IPW, which is p_ = ^ for the given values: like 
scattering amplitudes, the Wigner function sees particles 
of free mass to as having mass to, in an IPW. Away 
from the peak, the Wigner function becomes negative: 
as the shifted mass-shell condition cannot be solved, this 
naturally corresponds to quantum (off-shell) effects. For 
p_ > ip, there is in fact no real solution to the deformed 
mass-shell condition p 2 = M 2 (k.x, k.y) for any k.x or k.y. 



A. Monochromatic limit 

We now perform the y + integration in (|14[) and turn to 
the Wigner function itself, beginning with the IPW re- 
sult. To do so we add a Gaussian regulator exp[— e (k.y) 2 ] 



1 For simplicity, we rescalo to dimcnsionlcss variables in all our 
plots as follows: consider the exponent of (1141 . and then absorb 
the frequency k + into the co-ordinates and momenta, so we have 
k.x —> x + , k.y — s> y + , p + /k + — > p + , fc + p_/m 2 — > p-. With this, 
x + , y + and p± become dimcnsionless. 



FIG. 3. The Wigner function, [a.u.], plotted at p± = 0, 
p+ = 2, a§ = 25 as a function of p_. We clearly observe 
convergence for decreasing values of e (= 10~ 4 blue/dotted 
line, 10~ 5 purple/dot-dashed, and 10 -6 red/dashed) towards 
the final result (solid). The peak sits at the shifted mass-shell 
condition, and the photon emission and absorption numbers 
discussed in the text are also marked on the axes. 



However, for both sufficiently larger and smaller values 
of p_ , the Wigner function becomes positive again. The 
behaviour of the low p_ values in Fig. [3] is reminiscent of 
the edge and higher harmonic structure in Fig. [T] which 
results from multi-photon effects. Indeed, one may ver- 
ify that the edges at p_ = ip and in Fig. [3] are the 
solutions of the equation (reinserting k + and m factors) 

(p + nk) 2 =ml, (18) 

for n = 1 and n = 2. Hence, these edges correspond 
to resonances in which an off-shell particle can go on- 
shell, in the m* sense, by absorbing photons from the 
background. The jump at p_ = on the other hand, 
solves (fTS)) for n = — 1: it corresponds to a particle going 
on-shell by emitting a photon of laser momentum k^. 

It is possible to reveal the different processes contribut- 
ing to the Wigner function by Fourier decomposing the 
y + dependence in M 2 (as is done to expose the structures 
in Fig. [1] see @). One finds 



W A (x,p) 



4|p_ 



E 



dy + exp 



W \P+ 



4p_ 



x exp 



J„(^sinc(^B) ,(19) 



with J n Bessel functions of the first kind, and n gives 
the number of laser photons participating. A particularly 
neat illustration is shown in Fig. |4] Here, we use carte- 
sian momentum components, plotting at fixed po, as a 
function of longitudinal momentum p% , in order to give 
an alternative (and perhaps more familiar) view of our 
results, and a second set of parameters given in the cap- 
tion. With this choice of parameters, the shifted mass- 
shell condition p 2 = m 2 corresponds to p$ = ±| in the 
plots. The Wigner function is negative between these val- 
ues, and this is precisely where there is no real solution 
to the deformed mass-shell condition. 




FIG. 4. Top: Wa evaluated at p ± = 0, p = 5/4, a 2 = 1/2 
as a function of pz, [a.u.]. Quantum effects (negative values) 
occur when there is no solution to the deformed mass-shell 
condition. Bottom: Contributions from n = — 1 (emission), 
n = and n = 1 ... 3 (absorption) photon processes. 



The contributions from n = — 1 to n = 3 are shown 
in the bottom panel of Fig. 2J One sees that the n = 
term is responsible for the strong peaks at the effective 
mass condition, but is otherwise negative. The n > 
terms generate the edge structure: the locations of the 
first and second edges marked, as 'u' and 'v' in the plot, 
correspond to multi-photon conditions (|18l) for n = 1 and 
2. Hence the n of the Bessel expansion indeed gives the 
photon number. The n = — 1 process contributes largely 
to the quantum part. We turn now to pulse profiles. 



B. Pulse profiles 

We begin with some general properties of W4, which 
are determined entirely by M 2 . As seen in Fig. [21 a pulse 
profile brings in a dependence on k.x. Here, we focus on 
the most interesting case, where k.x is chosen to lie in 
the centre of our pulses. 

In a pulse, M rises rapidly over the first laser cy- 
cle, followed possibly by a brief period of oscillatory be- 
haviour and then a fall back to the rest mass squared. 
Consequently, we split the integral in (|14p into two pieces, 
corresponding to these different behaviours. It is natu- 
ral (but not necessary) to make this split at the edge of 
the pulse as seen by the y + integral, and so we label our 
two pieces "in" and "out" . (The change from "IPW-likc" 
to "fall-off" behaviour occurs deeper inside the pulse for 
tighter envelope functions.) The first part of the Wigner 
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function, W| n , is clearly always finite. If the pulse enve- 
lope varies slowly over most of the pulse duration then 
M will closely resemble the IPW effective mass; in this 
case, Wf 1 may be approximately written as with 
M 2 — > M 2 ilc and an infra-red cutoff in the y + integral at 
the edge of the pulse. W™ then typically exhibits rapid 
oscillations around the shape of the IPW result. 

The remaining part of the integral gives us W£ ut . By 
construction, M 2 — m 2 falls off as 1/k.y (plus subleading 
terms) in a pulsed background. This allows us to give an 
analytic expression for W£ ut : 

Wr = ci<5(p + - PT) + C2V X — , (20) 

p + - p°- s - 

where p° s ' is shorthand for the on-shell front energy 
po.s. _ ^2 _|_ m 2 ) /4p_. Special attention should be paid 
to the appearance of the free mass m in p° s \ The Cj 
are finite, pulse-profile dependent functions of momenta, 




FIG. 5. Wa in the pulse profile JT6J, at p = 5/4, p± = 0, a 2 = 
1/2, n = 8, 64 and 512 cycles (descending) with x + — nir, i.e. 
in the middle of the pulse. The free mass-shell condition is at 
p 3 = ±|, the IPW shifted shell is at p 3 = ±\. 



which we do not state explicitly here. Such a contribution 
is always found in W4 for pulsed fields as a consequence of 
the fields eventually turning off, or at least going to zero 
at large k.y. Consequently, we see that infra-red con- 
tributions to the Wigner function always yield the free 
mass-shell condition when we have a pulsed background 
field. Other mass scales may still play a role in determin- 
ing the detailed form of the Wigner function (as appears 
to be the case in scattering calculations [47j|). 

Let us illustrate all this with the pulse (fl"6"|). for which 

M 2 (nn,k.y)/ m 2 -l = a 2 ^^- (21) 

64 k.y 

when k.y > 27m, i.e. outside the pulse, confirming the 
falloff behaviour. The Wigner function is plotted in Fig. [5] 
for n = 8, 64 and 512 laser cycles. We choose all other 
parameters as in Fig. |4] to allow easy comparison. The 
most striking features are the peaks at the free mass-shell 
condition (p^ = ±f), independent of the pulse length. 
As the number of cycles increases, rapid oscillations de- 
velop around the free mass-shell, while two broad peaks 
form, and narrow toward the IPW shifted mass-shell at 
Pz = ±j. Between these broad peaks, a negative dip 
forms, which resembles that in Fig. |4] However, while 
this signals the reappearance of IPW behaviour for long 
pulses, the delta function peaks at the free mass-shell are 
always present. 

Hence, the shifted mass itself does not play quite the 
same role as it does in the IPW: for small numbers of cy- 
cles, the shifted mass-shell condition of the IPW is com- 
pletely absent in the pulsed W4. Even for long pulses, 
where broad peaks around the shifted mass-shell signify a 
"re-dressing" of the electron by the laser field, the Wigner 
function is always dominated by a delta function impos- 
ing the free mass-shell condition. 

Toward larger p% in Fig. [5J one sees a structure simi- 
lar to the multi-photon effects of the IPW results which 
exists independent of the pulse length. Because of the 

10 i 
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FIG. 6. The harmonic structure for n = 8 (dotted), 64 
(dash-dotted) and 512 (solid), compared with the IPW re- 
sult (dashed). IPW- like harmonic behaviour appears for a 
large number of cycles, though the free mass-shell is always 
present, seen here as the edge at p% — | at the very left of 
the plot. 
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FIG. 7. Wa in the 8-cycle pulse (|16|) . as a function of p_. 
Other parameters and ranges as in Fig. [3] to allow direct com- 
parison. The shifted mass-shell peak is completely absent. 



scales chosen, this is clearest in the n — 8 plot, however 
it is also present for other choices of n. In Fig. [5] we com- 
pare this structure in n = 8, 64 and 512 cycles with the 
IPW result. One clearly sees that this structure begins 
to closely resemble the IPW multi-photon structure for a 
large number of cycles: note, though, that at the left edge 
of Fig. [6j the free mass-shell delta function cuts through 
the IPW-like structure. This leads us to the conclusion 
that multi-photon effects remain in the game when we 
go to pulse profiles, and for very long pulses these learn 
about the shifted mass m„, as is clear from Fig. [51 al- 
though m continues to play a role. This is especially 
interesting in light of the analysis of the SLAC E144 ex- 
periment, which confirmed multi-photon effects, but did 
not unambiguo usly detect effects stemming directly from 
the mass shift [igf . 

We conclude this section with a final comparison be- 
tween the high intensity IPW result in Fig. O and the 
Wigncr function for a high intensity 8-cycle pulse, shown 
in Fig. [71 which further confirms the observations above. 
We display the same domain and range in both plots for 
easy comparison. While the large momentum behaviour 
is reminiscent of the IPW results, including a sign change 
near the first emission point, the shifted mass-shell at 
p_ = is completely absent. (The free mass-shell is of 
course present, but lies outside the plotted range.) 



C. Equal front time 

We consider now the equal lightfront time Wigner func- 
tion. To this end, we perform the p + integration in (|14[) . 
which corresponds to setting y + = 0. As the effective 
mass M deviates from the free mass m due to correla- 
tions in front time, one sees that all effective mass effects 
must vanish when y + = and M = m. Indeed, one finds 
that, independent of the pulse profile, W£ = l/4|p_| just 
as in the free theory. Hence, calculations of, say, num- 
ber density at different front times will always agree, and 
so the triviality of W£ seems to reflect the well-known 



statement that the plane wave cannot produce pairs. It 
is tempting to conjecture a relation between this result 
and the 'triviality' of the lightfront vacuum [33^ . 



D. Dirac spinors 

As we have seen, one benefit of the lightfront approach 
is that the electron-positron Wigner function is particu- 
larly simple. It may be the case that one would also 
like to consider the Wigner function for commutators of 
the constrained field ip(-), or of the constrained and dy- 
namical fields. All these interdependent Wigner func- 
tions can be neatly represented by the Wigner function 
defined from the commutator of the full Dirac spinor 
( | [ipi v \ I ) ■ Writing $4 for the integrand of the scalar 
t/4, see ([13)) . the Dirac field density is 



$ 1 



(2tt) 3 
x {4 + m) 



1 



2k.q 



1 



2k.q 



2 1 



(22) 



where we recognise the spin structure of the fermionic 
Volkov propagator. The Wigner transform is then per- 
formed just as in ([2j) and all commutators may be iden- 
tified by projection onto the relevant (+) and (— ) com- 
ponents. 



V. TWO APPROACHES TO PAIR 
PRODUCTION 

We briefly outline how the earlier results may be ex- 
tended to include the effects of vacuum pair production. 
To do so, we add a second plane wave, i.e. a function of 
k'.x where k'.k 7^ and vector structure as in (| 12|) . The 
collision of two such pulsed plane waves allows us to study 
vacuum pair production with both electric and magnetic 
fields, and with finite size effects in two directions, since 
pairs are only produced where the fields overlap. 

We consider two experimental setups. The first is the 
collision of two high intensity lasers, as envisaged for ELI 
Q . In this situation, a perturbative field expansion is not 
permissible. However, if we take the laser pulses to be 
almost parallel (pair production in parallel pulses is con- 
sidered in [HI), then k.k! « because the vector product 
obeys k.k' = 1 + e « 1, which gives us a small expan- 
sion parameter e. Alternatively, for two optical lasers, we 
have k.k' /m 2 ~ 10~ 12 independent of their collision an- 
gle, which gives us another small parameter with which to 
construct the Wigner function. Note that both of these 
approximations arc independent of the field strengths: 
hence the corresponding expansions are nonperturbative 
in the field intensities, allowing them to be used in the 
analysis of, for example, multiple ELI-strcngth lasers. 
The detailed analysis of these expansions will be a part 
of our future research. 
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The second setup is the collision of an intense optical 
laser with a high frequency, low intensity X-FEL beam. 
If the amplitude of the latter is sufficiently small, it may 
be treated in perturbation theory. To illustrate, we begin 
with the differential equation denning the Wigner func- 
tion (13: 



(§# + I4-to)W = 0, 
where the gauge invariant operators D and II are 



(23) 



Dp = - e J ds F^x - isd p )dp , 

-1/2 
1/2 

n M = Vn - ie J ds s F^ v (x - isd p )dp . 



(24) 



-1/2 



One now perturbs the field strength, — > F^ v + SF^, 
where, for us, F^ v represents the intense (optical) plane 
wave, and SF^ represents the weak X-FEL field. The 
above operators are then expanded as — > + 
etc. Writing W — > W + SW, one immediately obtains the 
equation 



(|# + J/L- m)5W = -{{5$ + 6J/L)W 



(25) 



This equation is solved by inverting the operator on the 
left hand side to give a Green's function. This may be 
shown to be a gauge invariant version of the well known 
Volkov propagator. Thus, the perturbation SW will de- 
scribe the emission of electron-positron pairs, dressed 
by the strong background, from the weak background 
source. Again, the detailed investigation of this result 
will be delayed for future work: we mention briefly that 
the first correction to the equal lightfront time Wigner 
function, <5Wg , is no longer trivial, which, from our pre- 
vious discussion, we take as a confirmation of pair pro- 
duction. The Winger function can be seen to depend on 
both k' .x (which is, say, x~ if we consider a head-on col- 
lision) and also the original lightfront direction k.x. In 
particular, the dependence on the strong field remains 
nonperturbative, as seen in the plane wave W4. Thus, 
this approach may also be used to study the role of the 
effective mass in vacuum pair production. 



VI. CONCLUSIONS 

We have constructed the Wigner function for both 
scalar and spinor particles in a background plane wave 



field exactly. This gives a precise definition of the effec- 
tive mass which extends the infinite plane wave mass- 
shift to arbitrary pulse profiles. The Wigner function 
itself clearly exhibits, in the infinite plane wave case, the 



shifted mass-shell condition p 



multi-photon ab- 



sorption/emission, and quantum effects where the effec- 
tive mass-shell condition has no solution, and the Wigner 
function becomes negative. 

Going to a pulsed profile, we have observed multi- 
photon structures in both long and short pulses. We 
have shown that the Wigner function in a pulse always 
contains, explicitly, a delta function imposing the free 
mass-shell condition, independent of the pulse length or 
geometry. This is in contrast to the infinite plane wave 
case, where the full, nonperturbative Wigner function 
knows only about the shifted mass-shell. 

Finally, we gave various approximations with which to 
calculate the Wigner function in a more complex field 
geometry describing colliding laser pulses. This will al- 
low us to study pair creation from the vacuum using the 
lightfront Wigner approach. 

There are many other extensions of the current work 
which one could consider. For example, transverse size 
effects are expected to become more important as one 
goes to higher intensities, since high intensity laser pulses 
are generated by tight focussing. 

Furthermore, recent calculations have shown that pair 
cascades, through highly focused laser interaction with 
single electrons, can take place at intensities as low as 
10 24 W/cm 2 [50, [51(. Consequently, if large numbers of 
photons are being taken from the beam, giving rise to sig- 
nificant beam depiction, it will become necessary to go 
beyond the background field approximation in describing 
the laser. Although the principles behind the inclusion of 
effects stemming from beam depletion are reasonably well 
understood, going beyond the background field approxi- 
mation still constitutes a major computational challenge. 



ACKNOWLEDGMENTS 



We are grateful to Rcinhard Alkofer and Thomas 
Heinzl for useful discussions. F. H. is supported by the 
DOC program of the Austrian Academy of Sciences, and 
by the FWF doctoral program DK-W1203. A. I. and 
M. M arc supported by the European Research Coun- 
cil under Contract No. 204059-QPQV, and the Swedish 
Research Council under Contract No. 2007-4422. 



[1] A. Ringwald, Phys. Lett. B 510 (2001) 107. 

[2] G. V. Dunne, Eur. Phys. J. D 55 (2009) 327. 

[3] T. Heinzl and A. Ilderton, Eur. Phys. J. D 55 (2009) 359. 



[4] M. Marklund and P. K. Shukla, Rev. Mod. Phys. 78 
(2006) 591. 

[5] T. Heinzl et al, Opt. Commun. 267 (2006) 318. 



9 



T. Heinzl and O. Schroeder, J. Phys. A 39 (2006) 11623. [31 
H. Gies, Eur. Phys. J. D 55 (2009) 311. [32 
S. S. Bulanov et al, Phys. Rev. Lett. 104 (2010) 220404. 
C. Harvey, T. Heinzl and A. Ilderton, Phys. Rev. A 79 [33 

(2009) 063407. [34 

G. V. Dunne and T. Hall, Phys. Rev. D 58 (1998) 105022. [35 

E. Brezin and C. Itzykson, Phys. Rev. D 2 (1970) 1191. 

S. P. Kim and D. N. Page, Phys. Rev. D 65 (2002) [36 
105002. 

A. D. Piazza, Phys. Rev. D 70 (2004) 053013. [37 

C. K. Dumlu and G. V. Dunne, Phys. Rev. Lett. 104 

(2010) 250402. [38 

H. Gies and K. Klingmiiller, Phys. Rev. D 72 (2005) [39 
065001. 

R. Alkofer et al, Phys. Rev. Lett. 87 (2001) 193902. [40 

D. B. Blaschke et al., Phys. Rev. Lett. 96 (2006) 140402. 

F. Hebenstreit et al., Phys. Rev. Lett. 102 (2009) 150404. [41 
H.-W. Lee, Phys. Rept. 259, 147 (1995). [42 
H. T. Elze, M. Gyulassy and D. Vasak, Nucl. Phys. B 

276 (1986) 706. [43 

H. Haug and A-P. Jauho, Quantum Kinetics in Transport [44 
and Optics of Semiconductors, Springer Series in Solid- 
State Sciences vol 123 (2007) (Berlin: Springer). [45 

I. Bialynicki-Birula, P. Gornicki and J. Rafelski, Phys. [46 
Rev. D 44 (1991) 1825. 

F. Hebenstreit, R. Alkofer and H. Gies, Phys. Rev. D 82 [47 
(2010) 105026. 

P. Levai and V. Skokov, Phys. Rev. D 82 (2010) 074014. [ 
C. K. Dumlu, Phys. Rev. D 79 (2009) 065027. 
A. M. Fedotov et al, arXiv: 1008.2098 [hep-ph]. 
P. Zhuang and U. W. Heinz, Annals Phys. 245 (1996) [49 
311. [50 
Ch. Kurtsiefer, T. Pfau and J. Mlynek, Nature 386 
(1997) 150. 

G. Breitenbach, S. Schiller and J. Mlynek, Nature 387 [51] 
(1997) 471. 

P. Bertet et al, Phys. Rev. Lett. 89 (2002) 200402. 



N. Sengupta, Bull. Math. Soc. (Calcutta), 44, 175 (1952). 
S. J. Brodsky, H. C. Pauli and S. S. Pinsky, Phys. Rept. 
301 (1998) 299. 

T. Heinzl, Lect. Notes Phys. 572 (2001) 55. 

T. Heinzl, A. Ilderton, Opt. Commun. 282 (2009) 1879. 

A. I. Nikishov and V. I. Ritus, Sov. Phys. JETP 19 (1964) 

1191 [Zh. Eksp. Teor. Fiz. 46 (1964) 1768]. 

A. I. Nikishov and V. I. Ritus, Sov. Phys. JETP 19 (1964) 

529. 

N. B. Narozhnyi, A. Nikishov, and V. Ritus, Zh. Eksp. 
Teor. Fiz. 47, 930 (1964). 

C. Bamber et al, Phys. Rev. D 60 (1999) 092004. 

L. S. Brown and T. W. B. Kibble, Phys. Rev. 133 (1964) 
A705. 

T. Heinzl, D. Seipt and B. Kampfer, Phys. Rev. A 81 
(2010) 022125 

D. Seipt and B. Kampfer, arXiv: 1010.3301 [hep-ph]. 

F. Mackenroth and A. Di Piazza, arXiv: 1010.6251 [hep- 
ph]. 

D. Volkov, Z. Phys. 94, 250 (1935). 

T. W. B. Kibble, A. Salam and J. A. Strathdee, Nucl. 
Phys. B 96 (1975) 255. 

M. Bonitz et al, Contrib. Plasma Phys. 39 (1999) 329. 

F. Mackenroth, A. Di Piazza and C. H. Keitel, Phys. 
Rev. Lett. 105 (2010) 063903. 

T. Heinzl, A. Ilderton and M. Marklund, Phys. Lett. B 
692 (2010) 250. 

K. T. McDonald, Prepared for Conference on Probing Lu- 
minous and Dark Matter: Adrian Fest, Rochester, New 
York, 24-25 Sep 1999 

A. Di Piazza et al, Phys. Rev. Lett. 103 (2009) 170403. 

G. Korn and P. Antici (Eds.), ELI Grand Chal- 
lenges Meeting Report, http://www.extreme-light- 
infrastructure.eu/reports.php (2009). 

N. V. Elkina et al. arXiv: 1010.4528 [hep-ph]. 



